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Abstract 

The problem of separation of uranium isotopes in a concurrent gas centrifuge 
is solved analytically. Separative power of the optimized concurrent gas cen¬ 
trifuges equals to 5U = 12.7(V/700 m/s)^(300 K/T)L, kg SWU/yr, where 
L and V are the length and linear velocity of the rotor of the gas centrifuge, 
T is the temperature. This formula well agrees with an empirical separative 
power of counter current gas centrifuges. The optimal value of the separative 
power is not unique on the plane {pw,Vz), where pw is pressure at the wall of 
the rotor and Vz is axial velocity of the gas. This value is constant on a line 
dehned by the equation PwVz = constant. Equations dehning the mass flux 
and the electric power necessary to support the rotation of the gas centrifuge 
are obtained. 

Keywords: gas centrifuge, isotope separation, diffusion in strong 

centrifugal field, separative power 
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1. Introduction 

Separation of isotopes has traditionally been carried out using 
gaseous centrifugation. Only recently the problem of sedimenta¬ 
tion of isotopes in liquids and solids due to diffusion of the isotopes 
in strong centrifugal fields attracts attention of specialists in con¬ 
nection with design of the centrifuges which can provide the cen¬ 
trifugal field ~ 10® g at the temperature achieving 500"C' |[Ij. The 
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diffusion process in the condensed matter differs essentially from 
the diffnsion in the gas. Theoretical gronnds of diffnsion in con¬ 
densed media are still nnder debates (see |2j and reference therein). 
Specific featnres of the sedimentation of liquid and solid mixtures 
in strong centrifngal fields and corresponding experiments were 
discussed in |31|T1,[5]. The attempts to interpret these experiments 
can be fonnd in |6| and references therein. In the present paper we 
discuss one of the unsolved key problem of the traditional process 
of separation of nraninm isotopes in the gas centrifnges (hereafter 
GC). 

The transport of gaseons isotopes in centrifngal fields of the or¬ 
der of 10® g is rather specific in compare with the transport of the 
isotopes in a condensed media. In so strong fields the density of 
the gas varies on ~ 10® — 10® orders of magnitnde per centimetre 
along radins. Convection in the GC plays crucial role because it 
is artificially excited to increase the separative power of the GC 
[7J. The relationship between the diffusive flux and the flux due 
to convection also varies on the same ~ 10® — 10® orders of magni¬ 
tude. This complicates the solution of the problem because it is 
necessary to solve all the system of Navier-Stokes equations and 
pins to solve the diffusion eqnation for the mixture of isotopes in 
extreme conditions of strong variation of density of the gas. On 
this reason the problem of the isotope separation is solved np to 
now only numerically. 

Separation of heavy isotopes in gas centrifuges is used for industrial pro¬ 
duction of enriched uranium starting with the middle of past century. In 
spite of old history, the problem of separative power of an optimized GC is 
not solved yet. Solution of this problem is interesting from practical point 
of view. Knowledge or estimate of the separative power of GC is important 
for development of new more efficient GC and important for experts dealing 
with the problem of proliferation prevention of the separation technology. 

Attempts to estimate the separative power of GC has been made firstly 
by Dirac [8]. He has shown that the separative power 5Umax of any GC can 
not exceed the value 


^ Urn.n.cr. 


iipDL 

2 ^ 2RT ' 


( 1 ) 


where pD is the density of uranium hexafluoride {UFq) times the coefficient 
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of self-diffusion for uranium isotopes . AM is the mass dif¬ 

ference between two uranium isotopes , R is the gas-law constant, T is gas 
temperature, L is the length of the rotor of GC, V is linear velocity of the 
rotor rotation. 

At the beginning of 1960’s an Onsager group from US developed a the¬ 
ory called the pancake approximation that reduced the problem of the GC 
gas dynamics to solution of linear elliptical equations of the sixth-order par¬ 
tial derivatives for two variables derivatives ISlllo]. This approach gave the 
following equation for the separative power of the GC 


6U = (0.038U - 11.5)L, kg ■ SWU/yr. (2) 

It is important to note that in contrast to eq. Q where the separative power 
increases as U^, in eq.([^ the separative power growths linearly with V. 

Experimental data collected with help of large amount of Russian GC 
shine new light on this question. According to m the separative power is 
dehned by the following empirical equation 


where L is measured in meters. Recently this result has been well conhrmed 
by more extended experimental data [T2] . 

In the empirical eq. ^ 6U ~ U^. This dramatically contradicts to simple 
theoretical arguments. At relatively small separation coefficient q dehned as 
a ratio of the concentration of in the product hux over the concentration 
in the waste hux the separative power equals to 

W = (4) 


where 6 = P/F is the ratio of the product mass hux P over feed mass hux F 
[7]. Radial separation in the centrifugal held gives the following dependence 
of g on U 


q = exp 


2RT 


(5) 


which unambiguously gives that 6U ~ V^. This dependence takes place in 
0 but do not agree with the experiment. For many years this problem 
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remains a challenge for specialists. Recently a new equation defining the 
separative power of GC has been proposed in [13] 

where V is measured in m/s, L is rotor length in meters, - some numerical 
coefficient. This equation already correctly reproduces the empirical law 
(|^. Unfortunately eq. (|^ has been obtained exploring weakly controlled 
assumptions and details of calculations are unknown. Therefore, the question 
about dependence of the separative power of the GC on the parameters 
remains opened up to now. 

It is necessary to stress that we discuss here the separative power of GC 
optimized on all parameters which can be controlled by a designer. The sep¬ 
arative power is the function of a lot of parameters SU(y, L, T, a, ai, 02 ,...), 
where the series of parameters ai includes for example pressure at the wall of 
the rotor, F,6, variation of temperature along the rotor ST and many other. 
Optimization of the GC is reduced to a search of maximum of this function 
at the variation of all the parameters «*. Such a search is performed for 
every series of V, L, T, a. Therefore, the separative power of the optimized 
GC depends only on the limited set of the parameters V, L, T, a. Such a 
formulation of the problem carries additional difficulties in the solution of the 
problem because it is necessary not only to calculate the separative power of 
the GC but additionally to optimize (to find maximal value) in relation to 
all possible parameters at fixed V, L, a and T. 

In this paper we give the solution of the problem of the separative power 
of the optimized concurrent GC from primary equations. This type of GC 
has been considered firstly in [S], where the separative power of this type of 
GC has been estimated as 

AMV^ o 

5U = 0.166 ■2v^pDL(^^)^ (7) 

which is only 66% less the maximal possible separative power given by ([^ 
and well agree with the arguments above in favour dependence of 6U . 

In this work we give analytical equation for the separative power 
of an optimized concurrent gas centrifuge for the first time. In 
contrast to the previous results we show that the correct solution 
gives another functional dependence of the separative power on V 
and L close to the empirical one given by eq. (|^. 
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The paper is organized as follows. In the second section we present the 
scheme of the concurrent centrifuge, basic equations and assumptions. In 
sec. 3 the solution is described in details. In sec. 4 the optimized separative 
power is calculated and hnally we discuss the solution in last sec.5. 


2. Isotope separation in the concnrrent centrifnge 


2.1. Hydrodynamics 

The schematic view on the concurrent GC is shown in £g. Here we 
use a scheme intentionally simplihed for the numerical analysis. The working 
gas UFq is load into the GC with feed flux F from the top (inlet) with the 
velocity independent on the radius r. We assume that the pressure p 
and density p correspond to the full hydrodynamic equilibrium in the radial 
direction. They are 


p = p^exp(^^{{^f-1)^, ( 8 ) 

and 

p = p^exp|^^((^)2-l)^, (9) 

where p^ and are the pressure and density at the wall of the rotor, 
7 = 1.067 - adiabatic index of the working gas and C = ■>/ 'yRT/p is the 
sound velocity, while p = 0.352 kg/mol is the molecular weight of UF^. The 
output of the gas occurs at the bottom of the rotor (outlet). To simplify 
the hydrodynamics we assume that the product P and waste W fluxes are 
separated by a concentric tube with radius r*, which provides the specihed 
ratio of the product flux to the feed flux 9 = P/F. In this case the flow 
lines are the straight lines parallel to the axis. Here we neglect second order 
effects affecting the velocity of the gas due to the viscous stresses and heat 
conduction. Therefore, Vz is constant everywhere in the GG. 

F is connected with the parameters of the gas as follows 

r r o .7 271 p^Vza‘^C\ 71/2 

Jo ^ 

The product flux P is defined as 


P = 


27ip^Vza^C^ ^ 7 I /2 

pVz2'Krdr = - — -(exp(^((-) - 1)) - exp(-^ 


))• 

( 11 ) 
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Hereafter we neglect exp(—2^) because its value is close 10 — 10 for 

typical parameters of GC. Then, radius r* is defined by the equation 

= ( 12 ) 

2.2. Separative power 

Concentration c in the feed flux equals to natural Cq = 7.114 ■ 10“^. We 
assume that c does not depend on r in the feed flux, c increases near the 
rotational axis and reduces at the wall upon motion of the gas along the 
rotor. The separative power of the GC is dehned as follows 

5U = PG{cp) + WG{cw)-FG{cf), (13) 

where G(c) = (2c— 1) In is the separative potential introduces by Dirac 
and Peierls [8]. The generalization of this equation on the case of the nonuni¬ 
form flow through the feed inlet, product and waste outlets gives [7] 

SU = j G(c)pvdS. (14) 

In our case the integration over the surface S which covers all the boundaries 
of the working volume gives 

6U = 2tt f pVz{G{c) — G{co))rdr, (15) 

Jo 

where G{c) is taken at the outlet. 


3. Solution of the problem 


3.1. Basic equations and assumptions 

The equation dehning diffusion of the mixture of uranium isotopes in 
gaseous UFq is as follows 

doc d.Jh 

where the components Jk of the flux of UF^ with the light uranium isotope 
in the cylindrical system of coordinates are given by the equations na 


Oc 

Jz = pVzC- pD — 


( 17 ) 
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and 


(18) 


dc AM 

Jr = pVrC - pD{— + —C(1 - C)^), 

where Vr is the radial component of the gas velocity, c is the concentration 
of UFq with light isotope . In the case nnder the consideration Vr = 0, 
p shonld be taken from eq. ([^. The natnral concentration cq satisfies to the 
condition cq -C 1. Concentration does not change strongly in the concnrrent 
GC. Therefore the condition c -C 1 is valid everywhere and we can consider 
the term (1 — c) as a constant eqnal to (1 — cq). Product pD does not depend 
on pressure. Therefore we consider it as a constant. 

Substitution of the components of J into (16) gives the following equation 
for the steady state diffusion 


dpVzC n ^ 

dz ~ 


rdc 


(19) 


where 5 = ^(1-Co). 

Let us introduce new variables t = exp - 1))’ ^ and 

express concentration c as 


c = Cot (t, z). 


( 20 ) 


In these variables density p = p^t^. For Y we obtain the following equation 
in these variables 


d dY oPwVzd, 


5‘^Y 


0 . 


( 21 ) 


Here we neglect variation of (r/a)^. This is reasonable because the geomet¬ 
rical scale Ar on which density, pressure and t vary along radius is small 
compared with a. Indeed, 


Ar _ 

~V ~ 7F2- 


( 22 ) 


For UF^ typical C = 86 m/s, while V > 600 m/s. This means that ^ < 
2 ■ 10-2. 


Let us express pw in (21) trough F using 
form 

d dY 2 2F , G2 dY 

7rapD^^^~^ ^ ^ 


. After that eq. (|21|) takes a 

(23) 


dz‘^ 


d‘^Y 

2"^ ^ -52y = 0. 
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Here it is convenient to replace variable 5 on variable ^ according to the 
equation 




7^2 


(24) 


Then we obtain 


d BY . F 

t—t— --— 

dt dt -RapD 




de 


(25) 


This equation can be simplihed. Distribution of Y is as follows at the inlet 


y((,o) = t‘. 


(26) 


Maximal possible variation of Y corresponds to an inhnitely long rotor, where 
the concentration of light isotope goes to the nonuniform distribution dehned 
by equation 

Y{t,oo)=t-^. (27) 

The main contribution into separative work gives the concentration at radius 
r* separating the product and waste fluxes because the mass flux density in 
the product flux is maximal here. The variation AY in this case is as follows 


AY = 9 


-5 




(28) 
1 we 


where 9 = is the value of at the r*. Taking into account that 9^ ~ 
have an estimate 

Ay^2(51n(^). (29) 


^9' 


Now let us take into account that ~ where The term 


can be neglected in eq. (25) provided that it is much less the last term 
6^Y in this equation. This condition gives 

( 30 ) 

The last condition is fulhlled for reasonable parameters of the centrifuges. 
This means that in the concurrent centrifuge we can neglect the axial diffusion 
in compare with the radial one provided that condition (30) takes place. The 
basic equation for the concentration takes the following form in this case 

F BY 


7^2 


B BY 




TiapD B^ 


6^Y = 0 . 


(31) 









3.2. Boundary conditions 

Boundary condition at the inlet (feed flux) follows to distribution (26). 
At the outlet end of the rotor no boundary conditions should be specihed 
because we neglect the axial diffusion and the concentration is advected to 
this boundary from the computational domain. Boundary condition at the 
wall of the rotor corresponds to zero flux of the concentration 


dc d In p 


(32) 


Substitution of the pressure distribution (|^ into this equation gives the fol¬ 
lowing boundary condition for function Y 


71/2 


dY 




(33) 


At the wall of the rotor [t = 1) this gives the following boundary condition 

(34) 


dY 

t— + 5Y = t). 
ot 


We formally extend the variation of t from 1 at the wall of the rotor to 
0, because at the axis t is of the order 10“^^ — 10“^'^. We also neglect the 
variation of r/a assuming that this parameter equals to 1. Therefore, the 
boundary condition at f = 0 should be 


dY 




(35) 


3.3. Solution 


We can expand the solution of eq. (31) on functions exp(—fc^) in a sum 
of series as follows 

^(^>0 = y^exp(-A;i0Ai(0, (36) 


where functions Si satisfy to the equation 


d dSi , Fki 

_ t _ - 4- f_ — 

tdt dt ^irapD 


52 

- = 0 . 


(37) 


General solution of this equation can be expressed as a linear combination 
of Bessel functions Js{^t) and J^siKt) 


Si — E.Js{\it) BJ-si^it), 


(38) 
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where E and B are some constants and A* = y JapU ' define these con¬ 
stants let ns firstly consider the bonndary condition at t = 0. According to 
[T^ the fnnctions Js{Xit) and behave near the point t = 0 as 




1 


r(5 + l)^ 2 


j.sm = 




r(l-5)' 2 




4r(2 - 5) 


)• 


(39) 

Tims, one of the fnnctions is regnlar in this point, while another diverges. 
Snbstitntion of the fnnction into bonndary condition (35) gives the ex¬ 
pression 


25 


,X. 


i\2 


r(i + 5)^2 


(40) 


which does not eqnal to 0. Snbstitntion of into (35) gives the expression 


2X2-5^2-25 

22-5r(2-h)’ 


(41) 


which goes to 0 at t —)■ 0. Thns, Si = BJ_s{Xit). Bonndary condition at the 
wall of the rotor gives 

X^^^^J^ + 5J.sX) = 0. (42) 

It follows from the properties of Bessel fnnctions that ( |15j) 

+ SJ-sit) = (43) 

Therefore, the bonndary condition at the rotor wall gives ns the eqnation 
defining the eigenvalne of the problem 


Ji-sX) — 0 . 


(44) 


Fonr first eigenvalnes are Aq = 0; Ai = 3.82; A 2 = 7; A 3 = 10.15. 

The eigenfnnction at A = 0 can be defined directly from eq. 
follows 

Soit) = BoX. 


(37) 


as 


(45) 


Finally the solntion takes the form 


Y{t,0 = BoX + J2B^ew{-hOJ-5Xt). (46) 

i=l 
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Coefficients Bi in the expansion of the solution into the sum of the series (36) 
are dehned from the boundary condition at the inlet = 0. This gives 


t^ = 




(47) 


i=l 


Eigenfunctions J^siKt) are orthogonal to each other and to t~^ with weight 
function t. Therefore, multiplication of this equation on and integration 
over t gives the following value of Bq 


Bo = l-6. 

Coefficients Bi are dehned from the equations 


B,= 


Jo 


(48) 


(49) 


It can be shown that for the specihed boundary conditions at t = 0 and at 
t = 1 we have 


t J-s{Kt)tdt — ^2+(5 ’ 


and 


Jls{\t)tdt = -Jlsi^i)- 


(50) 


(51) 


4. Optimized separation power 


It follows from eqs. and (10) that 


27ipVzrdr = 2Ftdt. 


(52) 


Therefore the separative power (15) can be presented as 


5U = F (G'(c) - G{co))2tdt, 


(63) 


where integration is performed over outlet. It is convenient to express ki as 

napD 


h = 


(54) 
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In this case the solution can be presented as 


Let us introduce new variable x as follows 

F = TrapDXj^x- (56) 

Then the solution takes a form 

c(Lx) = ( 1 + ^ _'^^ t'^exp(-^)J_^(A^t)j , (57) 

and the separative power becomes 

5U = napDXl^x [ {G{c) — G{cQ))2tdt. (58) 

Jo 


Substitution of (24) into this equation at z 
for 6U 


6U = 71 {- 


pDXl 


71/2 


: L gives the following expression 
L<I)((5,co,x), (59) 


where function 4) is 


®(^,Co,x) = X 



(G(c) - G{cQ))2tdt. 


(60) 


This function depends only on 3 parameters h, cO and y. Dependence of 4) 
on 5 and cO can be neglected at 5 ^ 1 and cO 1. These parameters 
are specified by the properties of the working gas. The only variable 
parameter is y. This is the dimensionless feed flux. Optimization of the 
concurrent centrifuge is possible only on this flux. This has clear physical 
sense. At the hxed length of the centrifuge L and pressure a small feed fluxes 
result into a small velocity Vz- Slow advection results into maximal possible 
separation of the isotopes. In this conditions decrease of the feed flux results 
into decrease of 5U because of the first term in the right hand side of eq. (53). 
In the opposite case of large feed flux, velocity Vz can be so large that the 
separation is negligible in compare with the advection of the concentration 
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to the outlet. Further increase of the feed flux results only to the decrease 
of the separation. At very large feed fluxes we also have to expect decrease 
of 6U because of decrease of the integral in the right hand side of eq. (53). 


Therefore, 6U have to have a maximum at some feed flux which correspond 
to the optimal feed flux in the centrifuge. The dependence of $ on y is shown 
in £g. This function has a maximum equal to 0.9 at X = 1. 

To express the optimized separative power of the centrifuges in conven¬ 
tional Separation Work Units (SWU) per year it is necessary to multiply eq. 
(59) on 1 year in seconds and on a factor ||| of ratio of weight of metallic 


uranium over weight of the working gas. After that we obtain 


5U = 12.7 


V 


700 m/s 


300 K' 
T 


L, kg-SWU/year 


(61) 


At the calculation of the separative power we assumed conventional value for 
pB = 2.3 ■ lO-'^ Pa • s. 


5. Discussion of the results 

In this work we have solved analytically the problem of diffusion 
of the uranium isotope mixture in the prescribed uniform flow of 
the gas with an exponential profile of density along radius in the 
concurrent centrifuge. The separative power of the GC has been 
calculated and optimized on the basis of this solution. For the first 
time an analytical expression for the optimized separative power of 
the gas centrifuge is obtained practically from the first principles 
which agree with experimental data (compare eqs. (|^ and ( |6lD ) . 
Equation (61) shows that in contrast to the result by Cohen |H] the optimized 
separative power of the concurrent gas centrifuge is proportional to V^L. 
This dependence well reproduces the experimental data for Russian counter 
current centrifuges approximated by eq. (|^. Surprisingly, the numerical 
coefficient in this dependence coincides with the experimental one in the 
limits of uncertanties. The separative power of the concurrent centrifuge 
does not depend on the rotor diameter. This is the only difference between 
equations ([^ and (61). 

It is interesting to understand why the optimized separative power of the 
concurrent centrifuge depends on the velocity of the rotor rotation as and 
this does not contradict to eq. (|^. It follows from eq. (55) that maximal 
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value of g — 1 equals to 


9-l = U"‘-l = 2Un-, 


(62) 


where = 6 (see eq. (12)). If we kept the separation of the fluxes 6 constant 


upon increase of V the value g — 1 remains constant as well. According to 


eq. (12) this means that upon growth of V the radius r* of the concentric 


tube separating the fluxes increases, keeping the product 


- 1 


(63) 


constant. As the result all the growth of 6U with V is due to growth of 
the feed flux. This flux increases with V because the rate of the radial 
separation increases as V^. Therefore, it is possible to increase the mass 
flux proportional to to provide efficient separation. The same is valid 
for dependence of 517 on L. g and 9 do not depend on L. The feed flux 
increases proportionally to L because the working gas has more traveling 
time for radial separation. 

The optimized separative power does not depend on the pressure of the 
gas at the wall of the rotor. This is another important feature of the obtained 
solution. It is clear why this occurs. Everywhere density enters into the 
equations in combination pvz- Product pD depends only on temperature. 


not density. Correspondingly, solution (55) depends on the feed flux F. 


Due to eq. (10) the optimal separative power is not unique on the plane of 


parameters Pw,Vz- Combining eqs. (10) and (56) at x = 1 we obtain that 


the optimal 6U is constant along the line located on the plane p^, defined 
by the equation 

IpDXlLV^ 


Pv 


4' 


C2 


(64) 


In the case of specification of pressure the axial velocity will be 


Vz = 60 


100 mm Hg^ f 6 cm 
Pw 


17 


700 m/s 


300 K 


L, m/s. (65) 


Optimal feed flux equals to 


F = 7rpDXl^L = 35 


17 


700 m/s 


L, g/s. 


( 66 ) 
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This means that the optimal concurrent centrifuge needs very large feed flux 
and provides very small coefficient of separation dehned by eq. (62). The 
electric power necessary to provide rotation of this feed flux can be estimated 
as 


1/2 

W = F— = 8.6 


V 


700 m/s 


300K 


L, kWatt. 


(67) 


Actually, this is the minimal necessary electric power spent for the work of 
the centrifuge. If to compare with the concurrent centrifuge, the production 
of the same separative work the concurrent centrifuge demands much more 
energy. This is the reason why the concurrent centrifuges are not explored 
for the industrial production of uranium isotopes. 
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Feed flux (F) 



Waste flux (W) 


Figure 1: A scheme of the concurrent centrifuge. The working gas is feed from the top of 
the rotor (inlet) with the distribution of pressure and velocity corresponding to the rigid 
body rotation of the gas. The product and waste fluxes leave the rotor at the bottom 
(outlet) and are separated by a concentric cylinder with radius r*. 
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